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Abstract 


The 3 He(^ _ ,n /i ) 3 H weak capture is studied using correlated-hyperspherical- 
harmonics wave functions, obtained from realistic Hamiltonians consisting 
of the Argonne or Argonne uis two-nucleon, and Tucson-Melbourne or 
Urbana-IX three-nucleon interactions. The nuclear weak charge and current 
operators have vector and axial-vector components with one- and two-body 
contributions. The axial-vector current includes the nucleon and A induced 
pseudo-scalar terms, with coupling constants gps and g* PS derived from pion- 
pole dominance and PCAC. The strength of the leading two-body operator is 
adjusted to reproduce the Gamow-Teller matrix element in tritium /3-decay. 
The calculated total capture rate is within ~0.5 % of the most recent mea¬ 
surement, 1496 ± 4 sec -1 . The predictions for the capture rate and angular 
correlation parameters A v , A t , and A a are found to be only very weakly de¬ 
pendent on the model input Hamiltonian. The variation of the observables 
with gps and g* PS and the theoretical uncertainties deriving from the model- 
dependent procedure used to constrain the axial current are investigated. 
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I. INTRODUCTION 


The /i weak capture on 3 He can occur through three different hadronic channels: 


t r + 3 He 
■- '- 3 He 
- 3 He 


/i 

/T 


3 H + z/ m (70%), 

n + d + t>fj (20%) , 

n + n +p + Up (10%) . 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 


The focus of the present work is on the first process. Some of the nuclear physics issues in 
muon capture have been reviewed recently in Ref. [1|. 

The reaction (|1.1| ) has been extensively studied through the years, both experimentally 
and theoretically. Measurements of the total capture rate have been performed since the 
early sixties U P t° uiitil recently. The latest very precise experimental determination 
of this observable ||, 1496 ±4 sec -1 , is consistent with the earlier measurements, the latter 
having considerably larger uncertainties, however. 

Theoretical studies of reaction HDD have been carried out within two different frame¬ 
works: the so-called “elementary particle method” (EPM) and the fully microscopic ap¬ 
proach. The EPM, first developed by Kim and Primakoff ||, is essentially a phenomeno¬ 
logical approach, which parameterizes the nuclear (charge-changing) weak current in terms 
of the trinucleon form factors, in analogy to the nucleon weak current, and then attempts 
to derive these from other experiments. Within the EPM, it was shown in Ref. |7] that, if 
the hyperfine structure of the p - 3 He system is taken into account and the direction of the 
recoiling triton can be detected, there are, in addition to the capture rate, other observables, 
i.e. angular correlation parameters, which are more sensitive than the capture rate itself to 
the value of the nucleon pseudo-scalar axial coupling constant gp$- Indeed, the possibility of 
determining g ps from measurements of muon capture observables is one of the motivations 
for the interest that this process has generated over the years. Recently, one attempt has 
been made to measure the angular correlation parameter A v [Q, though the experimental 
result, which to our knowledge represents the first significant measurement of this observ¬ 
able, is affected by large systematic uncertainties. Therefore, a comparison between theory 
and experiment will not be particularly meaningful for A v . Experimental results with an 
improved accuracy are highly desirable. 

In contrast, the fully microscopic approach is based on: i) 3 H and 3 He wave functions 
as accurate as possible, to reduce uncertainties related to nuclear structure; ii) a realistic 
model for the nuclear weak current and charge operators. The first microscopic calculation of 
reaction (El) was performed by Peterson in 1968 ||, and was reconsidered and improved by 
Phillips and collaborators in 1974 |]T0| . These studies, however, used nuclear wave functions 


which were approximate, and retained in the nuclear weak transition operators only single- 
nucleon terms, the impulse approximation (IA). 

In the early nineties, the muon capture on 3 He, including the total rate and angular cor¬ 
relation parameters mentioned above, have been extensively investigated, within the fully 
microscopic framework, by Congleton and Fearing 0, and Congleton and Truhlik ||12|| . 


The most significant improvements in these studies, relative to those of the late sixties || 
and early seventies [lTJ, are in the more accurate treatment of the trinucleon wave func¬ 
tions and of the weak interaction. The nuclear wave functions have been obtained from 
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a realistic Hamiltonian based on the Argonne vu two-nucleon 

three-nucleon [|lf|] interactions, using the rearrangement coupled-channel method |15 


|h| and Tucson-Melbourne 


The study in Ref. (TT| used the IA form of the nuclear weak current, and emphasized 
the need to go beyond single-nucleon contributions for a realistic description of the process. 
This next step was carried out in Ref. JT2j, where a model for two-body components in 
the nuclear weak current was explicitly constructed. The calculated capture rate is in good 
agreement with the measured value, although the theoretical prediction suffers from a 2 % 
uncertainty, which is rather large compared to the experimental error, and mostly arises 
from poor knowledge of some of the coupling constants and cutoff parameters entering the 
axial current. 

The present work sharpens and updates that of Ref. (T2j. Improvements in the modeling 
of two- and three-nucleon interactions and the nuclear weak current make the re-examination 
of process ( |1.1|) especially timely. The initial and final state wave functions have been ob¬ 
tained, using the correlated-hyperspherical-harmonics method, from a nuclear Hamiltonian 
which consists of the Argonne Uig two-nucleon [|l 6 j] and Urbana-IX three-nucleon [[17] inter¬ 


actions. To make contact with the study of Ref. |T2j , however, and to have some estimate 
of the model dependence of the results, the older Argonne U 14 two-nucleon and Tucson- 
Melbourne three-nucleon interaction models have also been used. Both these Hamiltonians 
reproduce the experimental binding energies and charge radii of the trinucleon systems. 

The model for the nuclear weak current used in the present work has been developed in 
Refs. |T3-|20|]. However, two additional contributions have been included: the one-body term 


associated with the induced pseudo-scalar charge operator of the nucleon, and the induced 
pseudo-scalar two-body term in the iVA-transition axial current. Both contributions are 
of order 0 (g 2 /m 2 ), where q is the momentum transfer in the process and m is the nucleon 
mass. They were neglected in the proton weak capture reactions studied in Refs. fT 8 |-pn 
for which q <C m. A brief description of these operators is given in Sec. HI. 


Some of the differences between the model for the nuclear weak current of Ref. IT 2 
and that adopted here should be noted. It is well known by now that the axial current 
associated with A excitation is the dominant (axial) two-body mechanism. I 11 the present 
work, its strength, i.e. the AA-transition axial coupling constant g \, has been determined by 
fitting the measured Gamow-Teller matrix element in tritium /3-decay. The inherent model 
dependence of this procedure has been shown to be very weak in studies of the proton 
weak captures on 1 H [|l^] and 3 He |^0[- I 11 Ref. |18] predictions for the 1 H(p,e + z/ e ) 2 H cross 


section, obtained with a variety of modern high-quality two-nucleon interactions, differed by 
significantly less than 1 %, once the coupling constant g\ had been fixed as described above 
within each given model Hamiltonian (for further discussion of this point as well as of the 
reasons for such a weak model dependence, see Ref. |]T 8 |). In Ref. | 12 |, on the other hand, 


g \4 is related to the nNA coupling constant and values ranging from the quark-model 

to the Skyrme-soliton model predictions are used for /Ava- 

There are additional differences in the detailed form of the pion range operators, which in 

lO 


Ref. 


were derived from a phenomenological chiral Lagrangian containing contributions 


from 7 r- and Ai-pole mediated currents. Congleton and Truhlik, though, ignored p-meson 
contributions both in the axial-vector and vector sectors of the weak current. These are re¬ 
tained in the present work. However, as it is clear from Ref. |12j and also from Sec. [TV| below, 


these differences have little numerical impact on the calculated muon capture observables. 
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Finally, the induced pseudo-scalar term in the A axial current is ignored in Ref. [12 


while here it is determined using pion-pole dominance and the partially-conserved-axial- 
current (PCAC) hypothesis. The induced pseudo-scalar coupling constant g* PS is related to 
g* A via the (extended) Goldberger-Treiman relation j2l] . 


A crucial issue, though, remains to be addressed, namely the extent to which the present 
model for the nuclear weak current is successful in predicting observed weak transitions 
(note that the cross sections of the proton weak capture processes mentioned above are not 
known experimentally). The present work fulfils this need by showing that the calculated 
rate for capture on 3 He is in excellent agreement with the measured value. 

This manuscript falls into five sections. In Sec. |I1| explicit expressions for the rate and 
angular correlation parameters are derived in terms of reduced matrix elements of multipole 
operators, while in Sec. H the model for the weak current is succintly described. The results 
are presented and discussed in Sec. |TV], and some concluding remarks are given in Sec. 0. 


II. OBSERVABLES 


The muon capture on 3 He is induced by the weak interaction Hamiltonian [22,23 


G 


H W =^J ,ix /..(X i/Vx i , 


( 2 , 1 ) 


where Gy is the Fermi coupling constant, Gy= 1.14939 xlO 5 GeV 2 24], and l a and j a are 
the leptonic and hadronic current densities, respectively. The former is given by 


L(x) = e 


—ik„x t— 


u(K, K) (1 - 7 5 )Vv(x, Su) , 


( 2 , 2 ) 


where ^(x, s M ) is the ground-state wave function of the muon in the Coulomb held of the 
3 He nucleus, and u(k u , h v ) is the spinor of a muon neutrino with momentum k*,, energy 
E u ( —k v ), and hclicity h u . While in principle the relativistic solution of the Dirac equation 
could be used, in practice it suffices to approximate 


Vv( x , Sp) ~ il>u{x)x{sn) = 'ipis(x)u(k fl , sj -*• 0 , (2.3) 

since the muon velocity v M ~ Za C 1 (a is the fine-structure constant and Z= 2). Here 
-01 s (x) is the Is solution of the Schrodinger equation and, since the muon is essentially at 
rest, it is justified to replace the two-component spin state y(s M ) with the four-component 
spinor ^(k^, s M ) in the limit k^ —> 0. This will allow us to use standard techniques to carry 
out the spin sum over s^ at a later stage. 

In order to account for the hyperkne structure in the initial system, the muon and 3 He 
spins are coupled to states with total spin / equal to 0 or 1. The transition amplitude can 
then be conveniently written as 

■'» (/, fz, 4. K) = ( 3 H, 4; v, h„ I H w I (^, 3 He); /, f z ) 

- E(Wbsl/.4)i4A„ sj ( 3 H,4b'4q)| s He, Ss ) , (2.4) 

S m ,S3 ^ ^ 


where 
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sfi) — ^(k^, h u ) 'jfj- (1 75)w(k At , s^) , 

and the Fourier transform of the nuclear weak current has been introduced as 

f( q)= / dxe iqx j CT (x) = (p(q),j(q)) , 


(2.5) 


( 2 . 6 ) 


with the leptonic momentum transfer q defined as q = — kj, ~ — k^. The Bohr radius 

of the muonic atom in the ground state is about 130 fin, i.e. much larger than the nuclear 
radius, and it is therefore well justified to factor out ^i s {x) from the matrix element of j CT (q) 


between the trinucleon ground states, by approximating it as [22,23 


Kr = Wn(o)r = n 


(2 am r 


7r 


(2.7) 


where r/q s (0) denotes the Bohr wave function evaluated at the origin for a point charge 2e, 
m r is the reduced mass of the 3 He system, and the factor 1Z approximately accounts 


for the finite extent of the nuclear charge distribution |22|,|23|. The value 1Z=0.98 is used 


here [BT 


Standard techniques [0 


are now used to carry out the multipole expansion of the 
weak charge (p(q)) and current (j(q)) operators in the general case in which 9 is the angle 
between the spin quantization axis (the z-axis) and the leptonic momentum transfer q: 


3 H,4|p(q)| 3 He,s 3 ) = V2tt Y V / 2/TTi 1 d l mC 

1 = 0,1 


-6)(^s 3 ,lm\~s' 3 )Ci(q) , (2.8) 


( 3 H,4|j 2 (q)| 3 He,s 3 ) = ~V2n Y V^T+li 1 d l mfi (-9) {^s 3 ,lm\^s 3 ) L t (q) , (2.9) 

i=o,i 1 A 


( 3 H,4|j A (q)| 3 He,s 3 ) = V3n i d^_ x (-0) (^s 3 , l + E^q)] , (2.10) 

where m=s' 3 — s 3 , A=±l, and Q, L;, Ei and Mi denote the reduced matrix elements 
(RME’s) of the Coulomb (C), longitudinal (L), transverse electric (E) and transverse mag¬ 


netic (M) multipole operators, as defined in Refs. [2(],^,^3|] . The d! mm , are rotation ma¬ 
trices in the standard notation of Ref. [^5[. Since the weak charge and current operators 
have scalar/polar-^vector (V) an d pseudo-scalar/axial-vector (A) components, each multipole 
consists of the sum of V and A terms, having opposite parity under space inversions m- 


Parity and angular-momentum selection rules restrict the contributing RME’s to C'o(V), 
C\(A), L 0 (V)i Li(A), E\(A) and Mi(V) in the 3 ffe(/i“, r / M ) 3 ff process. 

When the triton polarization is not detected, the differential capture rate for the reac¬ 
tion (|1.1|) is given by 


dT = 2n 8 ( m M + m T — E v — \Jnit + kl ) 


,—i2 dk^ 

(2tt) 3 


-w 


( 2 . 11 ) 


where m T1 and mt are the rest masses of the muon, 3 He, and 3 H, respectively, and the 
binding energy of the muonic atom has been neglected. Note that the following definition 
has been introduced: 
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|7V| 2 = E E p (/. A)I TwU,f*,s‘ z ,K 


( 2 . 12 ) 


s'o,h u fj z 


where P(f, f z ) is the probability of finding the /! 3 He system in the total-spin state / f z ). 
Integrating over the neutrino energy, the differential capture rate reduces to: 


^ i r 

d(cos6>) 2 0 


1 + A V P V cos 9 + A t P t cos 2 9 — — ^ + z4 A P A 


where the total capture rate T 0 reads 


r n = G\El (1-—) |^ av|2 

m t J 


Is I 0 ) 


with To denoting the following combination of R.ME’s 

To = \C 0 (V) - L 0 (V)\ 2 + I CM) - Li(A)\ 2 + I hhiV) - Ei(A)\ 2 
The angular correlation parameters A v , A t and 71 A are given by: 


A v — 1 + =- 


2 Im 


(C 0 (V) - L 0 (V))(C 1 (A) - L 1 (A))* - | M^V) - E 1 (A)\ 


(2.13) 


(2.14) 


(2.15) 


(2.16) 


A,= 


4 1 
3r; 
i 

'72 

1 


Im 

Im 


^=Re 

72 


2 1 


(C 0 (V) - LoiV^CM) ~ Li(A))* 
(C 0 (R) - L 0 (R))(M 1 (R) - E^A)) 
\Ci(A) - L l (A))(M l {V) - Ex(A)) 


-| Ml {V) - E 1 (A)\ 2 


(2.17) 


3 Tn 


72Im[(Co(R) - L 0 (V))(Mi(V) - E^A)) 
-T2Re[(C 1 (7l) - L 1 (A))(M 1 (V) - E X (A)) 
(C 0 (V) - LoiV^C^A) - L^A)) 


Im 


|Afi(U) - £' 1 (A)| L 


(2.18) 


Finally, the coefficients P v , Pt and P A are linear combinations of the probabilities P(f, /-), 
and are defined as |7|JTT11 


(2.19) 


P„ = P(1,1)-P(1,-1) , 

P t = P(l,l) + P(l,-l)-2P(l,0) , 

P a = P(1,1) + P(1,-1)+P(1,0)-3P(0,0) = 1 — 4 P(0, 0) 


Therefore, P v and Pt are proportional to the vector and tensor polarizations of the /=1 
state, respectively, while P A indicates the deviation of the /=0 population density from 
its statistical factor 1/4. Because of the small energy splitting between the /=0 and /=1 
hyperhne states (1.5 eV) compared to the jx~ 3 He binding energy, and hence small devia¬ 
tion of P(/, f z ) from its statistical value, direct measurements of the angular correlation 
parameters are rather difficult 


11 . 
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III. THE WEAK CHARGE AND CURRENT OPERATORS 


An exhaustive description of the model for the nuclear weak current has been recently 

H 


given in Ref. 


Here only its main features are summarized, and the new pseudo-scalar 


contributions are discussed. 

The nuclear weak current consists of vector and axial-vector parts, with corresponding 
one- and two-body components. The weak vector current is constructed from the isovector 
part of the electromagnetic current, in accordance with the conserved-vector-current (CVC) 
hypothesis. One important difference between the present calculations and those reported 
in Ref. |2(J is that the leptonic four-momentum transfer is not negligible, but in fact close 
to the muon rest-mass. Consequently, electromagnetic form factors need to be included in 
the expressions listed in Ref. The parameterization used for these reproduces available 
eN elastic scattering data. Furthermore, in the present work the Darwin-Foldy relativistic 
correction to the vector charge operator is also included. 

The one-body terms in the axial charge and current operators have the standard ex¬ 
pressions a obtained from the non-relativistic reduction of the covariant single-nucleon 
current, and include terms proportional to 1/m 2 , m being the nucleon mass. The induced 
pseudo-scalar contributions are retained both in the axial current and charge operators. In 
particular, the pseudo-scalar axial charge operator is taken as 


pIpsIt A ) — 


9ps 


2 mm. 


K - E ") 


■ q; a- 


(3.1) 


in the notation of Ref. [^U] . 

Again, because the leptonic momentum transfer involved in muon capture is not negli¬ 
gible, axial and induced pseudo-scalar form factors need to be included. These are parame¬ 
terized as 


9A{ql) 
9ps(ql ) 


9a 

(1 + qll^A? ’ 


2 m M m 

ml + ql 


9A{ql ) , 


(3.2) 

(3.3) 


where g 2 is the four-momentum transfer. The axial-vector coupling constant qa is taken to 
be [26] 1.2654±0.0042, by averaging values obtained from the beta asymmetry in the decay of 
polarized neutrons and the half-lives of the neutron and super-allowed 0 + —> 0 + transitions. 
The value for the cutoff mass A a is found to be approximately 1 GeV/c 2 from an analysis 
of pion electro-production data |27| and measurements of the reaction p(p M ,/r + )n |28[. The 
g 2 -dependence of gps is obtained in accordance with the partially-conserved-axial-current 
(PCAC) hypothesis, by assuming pion-pole dominance and the Goldberger-Treiman rela¬ 
tion |2l|^3|, m ^ here indicates the pion mass. 

Some of the two-body axial-current operators are derived from n- and p-meson exchanges 
and the p7r-transition mechanism. These mesonic operators, first obtained in a systematic 
way in Ref. |29|], give rather small contributions f2(J. The two-body weak axial-charge 


operator includes a pion-range term, which follows from soft-pion theorem and current al¬ 
gebra arguments 01 |3l|| , and short-range terms, associated with scalar- and vector-meson 


exchanges. The latter are obtained consistently with the two-nucleon interaction model, 


7 










































following a procedure 
current operators [^Oj] . 


[3^] similar to that used to derive the corresponding weak vector- 


The two-body axial charge operator due to iVA-transition is also 


included | 20 [, but its contribution is found to be very small. 


The dominant two-body axial current operator is that due to A-isobar excitation ||18|,|20 
We briefly review here its main features. The iVA-transition axial current is written as 
(notation as in Ref. ||20| ) 


J, ( 1 , (q;iV^A,.4) = - 


9l(d)S, + ■ q) 

2 mm,, 


3'qr.j 1 


i,± 


(3.4) 


where S* and T, are spin- and isospin-transition operators, which convert a nucleon into a 

has been obtained from a 

0 . 


A. The induced pseudo-scalar contribution, ignored in Ref. j20 


non-relativistic reduction of the covariant IVA-transition axial current 


The axial and pseudo-scalar form factors g A and g* PS are parameterized as 


9a(<£) = RAgA(ql) , 


dpsiVa) ~ 


2 m /t 771 

ml + ql 


aM 


(3.5) 


with gA^ql) given in Eq. (|T2|). The parameter R A is adjusted to reproduce the experimental 
value of the Garnow-Teller matrix element in tritium (3 decay, GT EXP = 0.957 ± 0.003 ||18|| , 
while the (^-dependence Q f g* ps i s again obtained by assuming pion-pole dominance and 


PCAC [0-11. The values for R A determined in the present study are listed in Table | for 
the four different combinations of interaction models. The experimental error on GT EXP is 
responsible for the 8-9 % uncertainty in R A . 

Before concluding this section, a couple of remarks are in order. First, it is important 
to note that the value of R A depends on how the A-isobar degrees of freedom are treated. 
In the present work, the two-body A-excitation operator is derived in the static A approxi¬ 
mation, using first-order perturbation theory (see Ref. |2(|]). This approach is considerably 
simpler than that adopted in Ref. 0 , where the A degrees of freedom were treated non- 
perturbatively, by retaining them explicitly in the nuclear wave functions |33]. The results 


for R 4 obtained within the perturbative (PT) and non-perturbative (TCO) schemes differ 
by more than a factor of 2-see Table VI of Ref. [2U|: i? J 4 (PT)= 1.22 and R j 4(TCO)=2.87Q. 
However, the results for the observables calculated consistently within the two different 
schemes are typically within 1 % of each other. 

Second, because of the procedure adopted to determine R A , the coupling constant 
g A =RA9A cannot be naively interpreted as the N A axial coupling constant. The excitation 
of additional resonances and their associated contributions will contaminate the value of g A . 
Indeed, the PCAC arguments used above imply g A / gA=fnNA/fnNN, where fnNN and f n NA 
are the 7 tNN and 7tN A coupling constants, and therefore one would obtain on the basis of 
Table | that (fnNA/fnNN) 2 is in the range 1.08-1.56, smaller than the value inferred from 
the A width, 4.67, and even smaller than the quark-model prediction, 2.88. 


1 Note that the value for ^^(PT) reported in Ref. [20] is slightly different from that given here in 
Table |[ since that value was obtained from a random walk consisting of 100,000 configurations, 
while the number of configurations sampled in the present work is 150,000. 







































IV. RESULTS 


In this section results for the 3 He(p - , z/ M ) 3 H capture process are reported. The trinucleon 
wave functions have been obtained from a realistic Hamiltonian consisting of the Argonne 
v 18 (AV18) [Jl6[| two-nucleon and Urbana IX (UIX) |l7j three-nucleon interactions. To com¬ 
pare with earlier predictions [Tl]JT^] for the same process, and to have some estimate of the 
model dependence, the older Argonne vu (AV14) [[13) two-nucleon and Tucson-Melbourne 
(TM) 0 three-nucleon interactions have also been used. Note that both the UIX and 
TM interactions have been adjusted to reproduce the triton binding energy. Finally, to 
investigate the effect of the three-nucleon interaction, predictions for muon-capture observ¬ 
ables have been made by including only two-nucleon interactions (AV14 or AV18) in the 
Hamiltonian models. 

The three-body bound-state problem has been solved with the correlated-hyperspherical- 
harmonics (CHH) method, as described in Refs. [33,35fl. It consists essentially in expanding 
the wave function on the CHH basis, and in determining variationally the expansion coeffi¬ 
cients by applying the Rayleigh-Ritz variational principle. 

The 3 H and 3 He binding energies are listed in Table [Qj for the different model Hamil¬ 
tonians employed in the present work. They are obtained including only the isospin 1/2 
components of the wave functions. These results, which are very accurate (the uncertainty 
is of the order of one keV), are in excellent agreement with the values calculated using other 
techniques (for a review, see Ref. [pp . 

Results for the capture rate T 0 and angular correlation parameters A v , A t , and A a, 
defined in Eqs. ( p,14|) - (|2.18|) , are presented in Table |ITJ. The uncertainty (in parethesis) in 
the predicted values is due to the uncertainty in the determination of the NA transition 
coupling constant g* A (see Sec. 0 and Table |) . The latter reflects the experimental error in 
the Garnow-Teller matrix element of tritium /1-decay. 

Inspection of Table 0 shows that the theoretical determination of the total capture rate 
To, when the AV18/UIX and AV14/TM Hamiltonian models are used, is within 1 % of 
the recent experimental result ||, 1496 ± 4 sec -1 . When the theoretical and experimental 
uncertainties are taken into consideration, the agreement between theory and experiment is 
excellent. Furthermore, the model dependence in the calculated observables is very weak: 
the AV18/UIX and AV14/TM results differ by less than 0.5 %. The agreement between 
theory and experiment and the weak model dependence mentioned above reflect, to a large 
extent, the fact that both the AV18/UIX and AV14/TM Hamiltonian models reproduce: 
i) the experimental binding energies as well as the charge and magnetic radii [[U]] of the 
trinucleons; ii) the Gamow-Tcller matrix element in tritium /3-decay. In this respect, it is 
interesting to note that the capture rates predicted by the AV18 and AV14 Hamiltonian 
models are about 4 % smaller than the experimental value, presumably because of the 
under-prediction of the binding energies and consequent over-prediction of the radii. This 
makes the relevant nuclear form factors, entering into the expression for the rate r 0 , smaller 
at the momentum transfer of interest, q ~ 103 MeV/c, than they would be otherwise. To 
study how the rate r 0 scales with the triton binding energy, we have repeated the calculation 
using CHH wave functions obtained with a modified AV14/TM Hamiltonian model, which 
gives for the 3 H and 3 He binding energies 9.042 and 8.349, respectively. The result for the 
rate T 0 is 1509 ± 7 sec -1 , while the angular correlation parameters are very close to the 
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AV14/TM values listed in Table |T|. Therefore, the rate To scales approximately linearly 
with the trinucleon binding energy. The values for the angular correlation parameter A v 
listed in Table 0 can be compared with the experimental result of Ref. |8|, 0.63 ± 0.09 
(stat.)lo!i 4 (syst.). Theory and experiment are in agreement, for any of the Hamiltonian 
models considered here. However, the experimental uncertainity is much larger than the 
theoretical one. 

The contributions of the different components of the weak current and charge operators 
to the observables and to the RME’s of the contributing multipoles are reported for the 
AV18/UIX model in Tables [TV], and [V| [VT|. respectively. The coupling constant g* A has been 
set equal to the central value of 1.17 (see Table |). The notation in Tables B0 and 0 is 
as follows. The column labeled “One-body no PS” lists the contributions associated with the 
one-body terms of the vector and axial charge and current operators, including relativistic 
corrections proportional to 1/m 2 . However, the induced pseudo-scalar contributions are 
not considered in both the axial current and charge operators. Therefore, the “One-body 
no PS”contribution is associated with the operators given in Eqs. (4.5)-(4.7), (4.8), (4.10), 
and (4.11)-(4.13) of Ref. [^], suitably modified by the inclusion of nucleon form factors, as 
explained in Sec. 0. The column labeled “One-body” lists the contribution obtained when, 
in addition, the induced pseudo-scalar axial charge and current operators, Eq. ( |3.1| ) and last 
term of Eq. (4.13) of Ref. H- respectively, are also included. 

The column labeled “Mesonic” lists the results obtained by including, in addition, the 
contributions from two-body vector and axial charge and current operators, associated with 
pion- and vector-meson-exchanges, i.e. the 7rV and pV for the vector current and charge 
operators, the 7 tA, pA and pnA for the axial current operator, and the 7 tA, sA and vA for 
the axial charge operator. We have used the notation of Ref. [^0|, where these terms are 
listed respectively in Eqs. (4.16)-(4.17), (4.30)-(4.31), (4.32)-(4.34) and (4.35)-(4.37). All 
these operators have been again modified by the inclusion of form factors. 

The column labeled “A no PS” lists the contributions arising from A excitation, but 
does not include those due to the induced pseudo-scalar A current of Eq. ( |3.4|) . The latter 
are retained in the column labeled “Full”. The associated operators are obtained, as men¬ 
tioned earlier in Sec. using perturbation theory and the static A approximation as in 
Eqs. (4.44), (4.48), (4.50) and (4.52) of Ref. [0. 


Note that in Tables [V] and the values for the RME Lq{V) have not been listed, 
since the charge and longitudinal multipole operators of the weak vector current, denoted 
respectively as Cn z (q\ V ) and L Uz {q\ R), are related via CYC as ['20 


Ui.W,V) = -- [H , Cu.(q-,V)] 


(4.1) 


In turn, this implies the following proportionality between the corresponding RME’s Co(R) 
and L 0 (V), L 0 (V ) = (m T — rn t — q 2 /2m t ) C 0 (V)/q, or L 0 (R) ~ —0.024 C 0 (R) for q ~ 103 
MeV/c. Finally, in Table [V] the induced pseudo-scalar axial contributions are present only 
in Ci(A) and Li(A), but not in E\(A), since the pseudo-scalar current is longitudinal. 

The importance of the induced pseudo-scalar contribution can be understood by inspec¬ 
tion of Table |TV]. The nucleon induced pseudo-scalar term in the axial current and charge 
operators reduce the value of To by about 16 %, while the changes in the polarization ob¬ 
servables are even larger. Far less important is the contribution from the pseudo-scalar 
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A-current, which reduces the value of To by less than 1 %. The changes in the polarization 
observables are also small, a few %. 

Among the observables, T 0 and A& are the most sensitive to two-body contributions in 
the weak current. These are in fact crucial for reproducing the experimental capture rate, 
see Table E3 Inspection of Table [V] shows that two-body contributions are significant in the 
RME’s Afi(y), h(A), and E 1 (A), but negligible in C 0 (V). The C 0 (V) and M X {V) RME’s 
are related by CVC to the corresponding RME’s of the isovector part of the electromagnetic 
current, since 


j-(q; v) 


r_,j iv (q;7) , 


(4.2) 


where j_(q; V) is charge-lowering weak vector current, ji v (q;7) is the isovector part of the 
electromagnetic current, and 77 is the (total) isospin-lowering operator. A similar rela¬ 
tion holds between the electromagnetic charge operator and its weak vector counterpart. 
Thus, if 3 He and 3 H were truly members of an isospin doublet, then the Cq(V) and M\[V) 
RME’s would just be proportional to the isovector combination of the trinucleon charge and 
magnetic form factors. Of course, electromagnetic terms and isospin-symmetry-breaking 
strong-interaction components in the nuclear potentials spoil this property. For example, 
the AV18/UIX model predicts for the isovector RME’s C' 0> i v (7) and Mi jiv ( 7 ) at q ~ 103 
MeV/c the values 0.3250 and -0.1385 (0.3254 and -0.1113 in impulse approximation), re¬ 
spectively. 

The C\(A) RME is about two orders of magnitude smaller than the leading RME’s, as 
expected on the basis of the following naive argument. The one-body axial charge density 
operator can be written approximately as (the notation is that of Ref. [j20[|) 


pl 1} (x; A) = ■ [pf, <5(x - rj) ] ~ (Ti ■ V;<5(x - r*) , 

2m 2m 


(4.3) 


where the term proportional to p; has been neglected, and the identity [A , B} + = [A , 7?]_ + 
2BA has been used. Here [ A , B\± denote the anticommutator (+) and commutator (-), 
respectively. We have also neglected the induced pseudo-scalar contribution. The one-body 
axial current density (its leading term) is 


jS,NR( x ; A ) = 9a cri 5(x - r*) , (4.4) 

and insertion of the approximation (|4.3| ) and Eq. ( |4.4| ) into the expressions for the charge 
and longitudinal multipole operators leads to the following relation between the associated 
RME’s: C\(A) ~ —(q/2m) Li(A), which, for q ~ 103 MeV/c, gives C\(A) ~ — 0.055 Li(A), 
i.e. the correct sign and order of magnitude obtained in the calculation. 

Lastly, from inspection of Table 0 . it is interesting to note that the contribution it A from 
the pion-exchange axial charge operator, which would be expected to be dominant among 
the two-body contributions to C'i(A), is also negligible. In fact, the operator structure 
of the corresponding C\(A) multipole is such that it cannot connect the dominant S-wave 
components in the 3 He and 3 H wave functions, and the associated matrix element is therefore 
highly suppressed. Furthermore, the ttA, pA and pi tA contributions to Li(A) and E\ (A) 
are about two orders of magnitude smaller than the leading one-body term (see Table |TV| ), 
and the relative signs of these contributions are such that they essentially cancel out in the 
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total sum. This feature of the mesonic contributions to the axial current was already found 


in other low-energy weak processes Hl8| |2C 


In order to compare with the results of Ref. pTl], the capture rate and angular cor¬ 


relation parameters have been calculated with the CHH wave functions corresponding to 
the AV14/TM Hamiltonian, and with a model for the nuclear weak current including only 
one-body terms. The values for the coupling constants and form factors entering the expres- 

imu 


sions for the charge and current operators have been taken from Ref. 


The comparison 


between the present and earlier predictions is shown in Table [VI Ij : there is satisfactory 
agreement between the two calculations. The remaining 1-3 % differences can presumably 
be explained as follows: i) the nuclear wave functions have been obtained with an AV14/TM 
Hamiltonian model with slightly different cutoff parameters |38|]; ii) the weak one-body op¬ 
erators in Ref. [[□]] include some of the next-to-next-to-leading orders in the non-relativistic 
expansion of the covariant single-nucleon current, proportional to 1/m 3 , these are ignored 
in the present calculation; iii) the numerical evaluation of the required matrix elements is 
performed with different techniques. Here, Monte Carlo methods based on the Metropolis et 
al. algorithm [[hj have been used. Typically, the statistical error on the calculated capture 
rate is less than 0.05 %. 

The results listed in Table H- column labeled “AV14/TM”, are also in good agreement 
with those of Table IX of Ref. ]T2|], although the treatment of the short-range behavior 
of the two-body terms in the weak current as well as the values for the vector and axial 
form factors, coupling constants, etc. in Ref. JT2| are slightly different from those adopted 


in the present work. It is important to emphasize, though, that the present model for the 
weak current reproduces well the available experimental data: i) the isovector component 
of the electromagnetic current, which by CVC is related to the weak vector current, leads 
to predictions for the isovector combination of the charge and magnetic form factors of 3 He 
and 3 H in excellent agreement with the measured values a 


of ~ 3 fm 


-l. 


up to momentum transfer 
ii) the two-body axial current operators are constrained to reproduce the 
Gamow-Teller matrix element in tritium /1-decay. 

To test the sensitivity of all the muon capture observables to the induced pseudo-scalar 
form factors gps and g * PS , Eqs. ( |3.3| ) and (|3.5|), we have repeated the calculation using 
AV18/UIX CHH wave functions and several different values of gps and g* PS in terms of their 
PCAC predictions g P ^ AG and g P F s CAG ■ We have assumed 


Rps - -ISc = • (4-5) 

9ps 9ps 

The variation of each observable in terms of Rps is displayed in Fig. |l|. The angular 
correlation parameters, in particular A t and Aa? are more sensitive to changes in g P g and 
g PS than the total capture rate, as first pointed out in Ref. [0. A precise measurement of 
these polarization observables could therefore be useful to ascertain the extent to which the 
induced pseudo-scalar form factors deviate from their PCAC values. 

Finally, by enforcing perfect agreement between the experimental and theoretical values, 
taken with their uncertainties, for the total capture rate To, it is possible to obtain an 
estimate for the range of values allowed for Rps- The procedure adopted is the following: 
i) we have considered the AV18/UIX minimum and maximum value for To (see Table ED), 
obtained with Rps= 1 and Ra—1- 08 and 1.26, respectively (see Table |). ii) For these two 
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values of Ra , we have tuned Rps to find r 0 within the experimental range. Our result for 
Rps is then 


Rps — 0.94 ± 0.06 


(4.6) 


This 6 % uncertainty is smaller than that found in previous studies [|lI],|T^,[4(]]. This sub 


stantial reduction in uncertainty can be traced back to the procedure used to constrain the 
(model-dependent) two-body axial currents described in Sec. m In this respect, it is inter¬ 
esting to note that ignoring altogether the mesonic axial contributions associated with the 
7T-, p- and /wr-exchange operators, and again re-adjusting the N A axial coupling constant 
to reproduce the tritium Gamow-Tcller matrix element (in this case, g* A = 1.32(9)^ is re¬ 
quired) lead to the following predictions for the muon capture rate and angular correlation 
parameters: T 0 =1479(7) sec^ 1 , A v =0. 5346(8), A t = 0.3666(14), and Aa=- 0.0988(13). In 
this case, the extracted value for the ratio Rps is 0.91 ± 0.06, in excellent agreement with 
the value of Eq. (|4.6|), suggesting that Rps is not too sensitive to these mesonic contributions. 


V. SUMMARY AND CONCLUSIONS 


Muon capture observables for the process 3 He(/i _ , i/ m ) 3 H have been calculated with very 
accurate CHH wave functions corresponding to realistic Hamiltonians, the AV18/UIX and 
AV14/TM models, and with a nuclear weak current consisting of vector and axial-vector 
parts with one- and two-body terms. The conserved-vector-current hypothesis has been 
used to derive the weak vector charge and current operators from the isovector electromag¬ 
netic counterparts, while the axial current has been constructed to reproduce the measured 
Gamow-Tcller matrix element of 3 H /3-decay. The axial current also includes the nucleon 
and A induced pseudo-scalar current operators. It should be emphasized that the model 
adopted for the electromagnetic current provides an excellent description of the 3 He and 3 H 
charge and magnetic form factors [[57J at low and medium values of momentum transfers. 

The predicted total capture rate is in agreement with the experimental value, and has 
been found to have only a weak model dependence: the AV18/UIX and AV14/TM results 
differ by less than 0.5 %. The weak model dependence can be traced back to the fact 
that both Hamiltonians reproduce the binding energies, charge and magnetic radii of the 
trinucleons, and the Garnow-Teller matrix element in tritium /3-decay. 

It is important to note that, if the contributions associated with two-body terms in 
the axial current were to be neglected, the predicted capture rate would be 1316 (1318) 
sec -1 with AV18/UIX (AV14/TM), and so two-body mechanisms are crucial for reproduc¬ 
ing the experimental value. The present work demonstrates that the procedure adopted 
for constraining these two-body contributions leads to a consistent description of available 
experimental data on weak transitions in the three-body systems. It also corroborates the 
robustness of our recent predictions for the cross sections of the proton weak captures on 
X H 


and 3 He 


[L9],^0|, which were obtained with the same model for the nuclear weak 

current. 

Finally, it would be interesting to study the 3 He(/i _ ,z/ At )nd and 3 He(/i“,z/ At )?mp pro¬ 
cesses, both of which have been investigated experimentally in Ref. |42|] and theoretically in 
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Ref. . Since the CHH method is suitable to solve for the three-body bound and scatter¬ 


ing states 0 , the study of these two processes is also possible. Work along these lines is 


vigorously being pursued. 
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TABLES 


TABLE I. Values for Ra = g\/gA, where g* A is the N A transition axial coupling constant (see, 
however, Sec. Ill for a discussion of the proper interpretation of g A ). The results are obtained by 
reproducing the experimental value of the Gamow-Teller (GT) matrix element in tritium /3-decay 
with CHH wave functions corresponding to the AV18, AV14, AV18/UIX, and AV14/TM Hamilto¬ 
nian models. The theoretical uncertainties are due to the experimental error with which the GT 
matrix element is known. 


Interaction Model 

9a/9a 

AV18 

1.25T0.10 

AV14 

1.11T0.09 

AV18/UIX 

1.17T0.09 

AV14/TM 

1.04T0.09 


TABLE II. Binding energies in MeV of 3 He and 3 H calculated with the CHH method using 
the AV18, AV14, AV18/UIX, and AV14/TM Hamiltonian models. Also listed are the experimental 
values. 


Interaction Model 


3 He 


3 H 

AV18 


6.917 


7.617 

AV14 


7.032 


7.683 

AV18/UIX 


7.741 


8.473 

AV14/TM 


7.809 


8.485 

EXP 


7.72 


8.48 

TABLE III. Capture rate To in sec 1 

, and angular correlation parameters A v , 

A t , and A a, 

as defined in Eqs. (2. 14|) (|2.18|), calculated 

using CHH wave 

functions corresponding to the AV18, 

AV14, AV18/UIX, 

and AV14/TM Hamiltonian models. The theoretical uncertainties, shown in 

parenthesis, reflect the uncertainty in the determination of the N A transition axial coupling con- 

stant g* A . 





Observable 

AV18 

AV14 

AV18/UIX 

AV14/TM 

r 0 

1441(7) 

1444(7) 

1484(8) 

1486(8) 

A v 

0.5341(14) 

0.5339(14) 

0.5350(14) 

0.5336(14) 

A t 

-0.3642(9) 

-0.3643(9) 

-0.3650(9) 

-0.3659(9) 

Aa 

-0.1017(16) 

-0.1018(16) 

-0.1000(16) 

-0.1005(17) 


17 























TABLE IV. Cumulative contributions to the capture rate Tq (in sec -1 ) and angular correlation 
parameters A„, A t . and A&- The CHH wave functions are obtained using the AV18/UIX Hamil¬ 
tonian model. The column labeled “One-body-no PS” lists the contributions associated with the 
one-body vector and axial charge and current operators, but no induced pseudo-scalar axial term is 
included. This is done in the column labeled “One-body”, while the column labeled “Mesonic” lists 
the results obtained by including, in addition, the contributions from meson-exchange mechanisms. 
Finally the column labeled “A-no PS” lists the results obtained by including also the A-excitation 
contributions, with g* A /dA set to the central value of 1.17 (see Table H), but excluding the A 
pseudo-scalar term, which is included in the column labeled “Full”. 


Observable 

One-body no PS 

One-body 

Mesonic 

A no PS 

Full 

r 0 

1530 

1316 

1384 

1493 

1484 

Ay 

0.7735 

0.5749 

0.5511 

0.5438 

0.5350 

A t 

-0.0840 

-0.3565 

-0.3679 

-0.3525 

-0.3650 

Aa 

-0.1424 

-0.0686 

-0.0810 

-0.1038 

-0.1000 


TABLE V. Cumulative contributions to the reduced matrix elements (RME’s) Co(V), Ci(A), 
Li(A), Ei(A), and M\ {V). The CHH wave functions are obtained using the AV18/UIX Hamilto¬ 
nian model. Note that Cq{V) is purely real, while the other RME’s are purely imaginary. Notations 
as in Table IV. 

RME 

One-body no PS 

One-body Mesonic 

A no PS 

Full 

C 0 (V) 

Ci (A) 

0.3280 

-0.7532 x 10 -2 

0.3277 

-0.4076 x 10 -2 -0.4135X 10 -2 

-0.4397x 10 -2 


Li(A) 

0.4058 

0.2590 0.2618 

0.2804 

0.2737 

Ei (A) 

0.5519 

0.5563 

0.5813 


M x (V) 

-0.1128 

-0.1314 

-0.1355 



TABLE VI. Individual mesonic contributions to the reduced matrix elements (RME’s) C'o(H), 
Ci(A), Li(A), Ei(A), and M\(V). The CHH wave functions are obtained using the AV18/UIX 
Hamiltonian model. Note that Cq(V) is purely real, while the other RME’s are purely imaginary. 
Notations as explained in the text. 


RME 

tt(V/A) 

P(V/ A) 

pirA 

sA 

vA 

C 0 (V) 

-0.3285 x 10 -3 

-0.6950X 10 -4 




Ci (A) 

-0.3253X 10 -5 



-0.2730 x 10 -3 

0.2174x 10 -3 

Li(A) 

0.2324X 10 -2 

-0.2894x 10 -2 

0.3409X 10 -2 



Ei (A) 

0.2539X 10 -2 

-0.4208X 10 -2 

0.6056X 10 -2 



Mi(V) 

-0.1597X 10 -1 

-0.2627x 10 -2 
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TABLE VII. Capture rate T 0 (in sec -1 ) and angular correlation parameters A v , A tl and 
A/± obtained with AV14/TM CHH wave functions, and only one-body operators (column labeled 
“One-body”) are compared with the results of Table 3 of Ref [111. 


Observable 

One-body 

Ref. |11] 

r 0 

1287 

1304 

A v 

0.579 

0.568 

At 

-0.351 

-0.356 

Aa 

-0.070 

-0.076 
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FIG. 1. Variation of the capture rate To and angular correlation parameters A v , At, and A/± 
with the induced pseudo-scalar coupling gps . The AV18/UIX CHH wave functions are used. For 
each observable, the ratio between the result obtained with the given value of gps and the PCAC 


prediction, listed in Table III, is plotted versus the ratio gps/9 ps" i = 9*ps!9*PS AC ’ see text). 
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